Contact mechanics: contact area and interfacial separation from small contact to full 

contact 
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Institut fur Festkorperforschung, Forschungszentrum Jiilich, D- 52425 Jillich, Germany 

We present a molecular dynamics study of the contact between a rigid solid with a randomly 
rough surface and an elastic block with a flat surface. The numerical calculations mainly focus on 
the contact area and the interfacial separation from small contact (low load) to full contact (high 
load). For small load the contact area varies linearly with the load and the interfacial separation 
depends logarithmically on the load. For high load the contact area approaches the nominal contact 
area (i.e., complete contact), and the interfacial separation approaches zero. The numerical results 
have been compared with analytical theory and experimental results. They are in good agreement 
with each other. The present findings may be very important for soft solids, e.g., rubber, or for very 
smooth surfaces, where complete contact can be reached at moderate high loads without plastic 
deformation of the solids. 
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1. Introduction 

What happens at the atomic and molecular level when 
surfaces come into contact with each other? And how do 
these events relate to macroscopic properties and obser- 
vations? These questions, which center on the phenom- 
ena of adhesion and friction, pose challenges not only in 
engineering but also in many areas of physical and bio- 
logical sciences [1]. Finding correlations and models that 
connect the atomic and macroscopic worlds usually is 
not easy. However, recently we have surprisingly found 
that the pressure distribution obtained from molecular 
dynamics calculations is in a good agreement with the 
prediction based on continuum contact mechanics, and 
in particular with the analytical contact mechanics the- 
ory of Persson0,i,[3. 

When two elastic solids with rough surfaces are 
squeezed together, the solids will in general not make con- 
tact everywhere in the apparent contact area, but only 
at a distribution of asperity contact spots [5|, [y]. The sep- 
aration iz(x) between the surfaces will vary in a nearly 
random way with the lateral coordinate x = (x, y) in the 
apparent contact area. When the applied squeezing pres- 
sure increases, the contact area A will increase and the 
average surface separation u = (u(x)) will decrease, but 
in most situations it is not possible to squeeze the solids 
into perfect contact corresponding to u = 0. The area 
of real contact, and the space between two solids has a 
tremendous influence on many important processes. 

Most studies of contact mechanics have been focused 
on small load where the contact area depends linearly on 
the loadBi&HEl,!!!. However, for soft solids, such 
as rubber or gelatin, and for smooth surfaces nearly full 
contact may occur at the interface, and it is of great inter- 
est to study how the contact area, the interfacial surface 
separation and stress distribution vary with load from 
small load (where the contact area varies linearly with 



the load), to high load (where the contact is (nearly) com- 
plete) . Here we will present such a study using molecular 
dynamics, and we will compare the numerical results with 
the prediction of the analytical contact mechanics theory 
of Persson. Our multiscale molecular dynamics approach 
PJJ] has been developed to study contact mechanics for 
surfaces with roughness on many different length scales, 
e.g., self-afflne fractal surfaces. 

This paper presents an extension of the work reported 
in two short publications 0, [l3j]. In Sec. 2 we briefly 
review the contact mechanics of Persson. In Sec. 3 we 
consider the relation between interfacial separation and 
squeezing pressure. Sec. 4 deals with the molecular dy- 
namics (MD) model. In Sec. 5 and 6 we compare the 
numerical results of the MD model with the analytical 
theory for the real contact area and interfacial separa- 
tion, respectively. In Sec. 7 we compare Persson theory 
with finite element calculations. Sec. 8 deals with how 
average surface separation depends on the magnification 
£. In Sec. 9 we consider the adhesion between randomly 
rough surfaces. Sec. 10 contains the summary and con- 
clusion. 

2. Theory: Contact area 

We consider the frictionless contact between elastic 
solids with randomly rough surfaces. If z = /ii(x) and 
ft-2(x) describe the surface profiles, E\ and E2 are the 
Young's elastic moduli of the two solids and v\ and V2 
the corresponding Poisson ratios, then the elastic con- 
tact problem is equivalent to the contact between a 
rigid solid (substrate) with the roughness profile h(x) — 
hi{x) + /&2(x), in contact with an elastic solid (block) 
with a flat surface and with an Young's modulus E and 
Poisson ratio v chosen so that [lfl \v\ 
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Persson dm has developed a contact mechanics theory 
where the surfaces are studied at different magnification 
( = Aq/A, where Aq is some reference length, e.g., the 
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FIG. 1: An rubber block (dotted area) in adhesive contact 
with a hard rough substrate (dashed area). The substrate 
has roughness on many different length scales and the rubber 
makes partial contact with the substrate on all length scales. 
When a contact area is studied at low magnification it appears 
as if complete contact occur, but when the magnification is 
increased it is observed that in reality only partial contact 
occur. 



roll-off wavelength of the surface roughness power spec- 
tra (see below), and A the shortest wavelength roughness 
which can be observed at the magnification see Fig. [TJ 
We define go = 27r/Ao. In this theory [2[ the stress dis- 
tribution P(cr, at the interface between the block and 
the substrate has been shown to obey (approximately) 
a diffusion-like equation where time is replaced by mag- 
nification and spatial coordinate by the stress a. When 
the magnification is so small that no atomic structure 
can be detected, the surface roughness will be smooth 
(no abrupt or step-like changes in the height profile) and 
one can then show [3] that in the absence of adhesion 
P(0, C) = 0. Using this boundary condition the solution 
to the diffusion- like equation gives the pressure distribu- 
tion at the interface (a > 0): 

= ^Tli (e- (CT - p)2/4G - c-(^)V«) , (1) 
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where p is the nominal squeezing pressure, and where 
n ( E \ 2 i-Cio 
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dq q 3 C(q) . 
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where is the smallest surface roughness wave-vector 
which may be of order 27r/L, where L is the linear size of 
the system. The surface roughness power spectrum [19[ 
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where (..) stands for ensemble average. The relative con- 
tact area 
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da P(a,C)=P(q), 
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FIG. 2: An elastic block squeezed against a rigid rough sub- 
strate. The separation between the average plane of the sub- 
strate and the average plane of the lower surface of the block 
is denoted by u. Elastic energy is stored in the block in the 
vicinity of the asperity contact regions. 



where q = £<2o- Sometimes, when we want to emphasize 
that P(q) depends on the pressure p, we will denote it by 
P p (q). Note that there is a delta- function contribution to 
P(<7, of the form [(Ao — A)/Ao]S(a) which arises from 
the non-contact area. Including this delta-function the 
integral of P(a, () over all a will be unity as expected for 
a probability distribution. In what follows we will always 
consider a > in which case P(crX) is given by (1). 
Substituting (pQ) into ([3]) gives after some simplifications 

— = 1 f da e-° 2 '* G = erf (-?—) U) 

Thus, for small nominal squeezing pressure p <IC G 1 / 2 we 
get 



(5) 
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A critical discussion of the theory presented above was 
given by Manners and Greenwood [14( . See also Carbone 
and Bottiglione[15|. 

3. Theory: Interfacial surface separation 

The space between two solids has a tremendous influ- 
ence on many important processes, e.g., heat transfer [2o|. 
contact resistivit y 12 lL lubrication [22|, sealing [23] and op- 
tical interference [24J. One of us has recently presented a 
simple theory for the (average) separation wasa function 
of the squeezing pressure p [4|]. The theory shows that 
for randomly rough surfaces at low squeezing pressures 
p ex exp(— u/uo) where the reference length uo depends 
on the nature of the surface roughness but is independent 
of p, in good agreement with experiments^. 

Consider an elastic block with a fiat surface squeezed 
against a hard rough substrate surface, see fig. [2l The 



separation between the average surface plane of the block 
and the average surface plane of the substrate is denoted 
by u with u > 0. When the applied squeezing force p 
increases, the separation between the surfaces at the in- 
terface will decrease, and we can consider p = p(u) as a 
function of u. The elastic energy U e \{u) stored in the sub- 
strate asperity-elastic block contact regions must equal 
the work done by the external pressure p in displacing 
the lower surface of the block towards the substrate, i.e., 



du f A p(u f ) = U e \(u), 



or 



p(u) 



1 dU el 
An du 



(6) 



(7) 



where An is the nominal contact area. Equation ([7j) is 
exact, and shows that if the dependence of the surface 
separation u on the squeezing pressure p is known, e.g., 
from finite element calculations or molecular dynamics, 
one can obtain the elastic energy U e \ stored in the as- 
perity contact regions [10]. This is an important result as 
U e \(u) is relevant for many important applications. 

Theory shows that for low squeezing pressure, the area 
of real contact A varies linearly with the squeezing force 
pAn, and that the interfacial stress distribution, and the 
size-distribution of contact spots, are independent of the 
squeezing pressure 25[ [26| . That is, with increasing p 
existing contact areas grow and new contact areas form 
in such a way that in the thermodynamic limit (infinite- 
sized system) the quantities referred to above remain un- 
changed. It follows immediately that for small load the 
elastic energy stored in the asperity contact region will 
increase linearly with the load, i.e., U e \{u) = unAnp(u), 
where un is a characteristic length which depends on the 
surface roughness (see below) but is independent of the 
squeezing pressure p. Thus, for small pressures (7) takes 
the form 

dp 

p(u) -uo-tz, 
du 



or 



p{u) oc e 



-u/uo 



(8) 



in good agreement with experimental data for the con- 
tact between elastic solids when the adhesional interac- 
tion between the solids can be neglected 0. We note 
that the result (8) differs drastically from the predic- 
tion of the Bush et al theory 0, and that of Greenwood 
and Williamson theory (GW)[8], which for low squeez- 
ing pressures (for randomly rough surfaces with Gaussian 
height distribution) predict p(u) oc u~ a ex.p(—bu 2 ), where 
a = 1 in the Bush et al theory and a = 5/2 in the GW 
theory. Thus these theories do not correctly describe the 
interfacial spacing between contacting solids. This is not 



surprising because these approaches assume a rigid sub- 
strate surface covered with flexible asperities. In reality, 
the bulk of the solids whose surfaces are in contact is not 
rigid. Furthermore, there exist a hierarchy of asperities 
on many length scales, all of which can distort. 

The elastic energy U e \ has been studied in Ref. [3] and 
12611. Here we will use 



C/ e i« A E*- / dqq 2 W(q,p)C(q), (9) 

Z Jqo 

where we have chosen q$ = q^ and where E* = E / (1 — 
v 2 \ q\ is the largest surface roughness wave vector and 

W(q,p) = P p (q) [7 + (1 - 7)f»] , (10) 
where P p (q) is given by Eq. (j4|): 



P P (q) = erf 



2G 1 / 2 (() 



(11) 



For complete contact (infinite squeezing pressure) P p = 1 
and thus W(q,p) = 1 and in this limit ([9]) is exact. For 
small squeezing pressure W(q,p) ~ jP p . The parame- 
ter 7 is of order w 0.4 (see below), and takes into ac- 
count that the elastic energy stored in the contact region 
(per unit surface area) in general is less than the average 
elastic energy (per unit surface area) for perfect contact 
(see Ref. |26|). The particular way we interpolate be- 
tween the limits W(q,p) = 1 for complete contact and 
W(q,p) = jP p for very small contact using (|T0)) was de- 
signed to give good agreement between the calculated 
[using (7)] interfacial separation, and the interfacial sep- 
aration obtained from Molecular Dynamics (MD) and Fi- 
nite Element Method (FEM) (see Sec. 6 and 7). 
Let us write (11) as[2l l27j 



2 f s (^P 2 
P P (q) = ~r dx e~ x , 
v n Jo 

where s(q) = w(q)/E* with 

w{q)=L f dq' q' 3 C(q') 
Using (12) gives 

du \J 7T du 



(12) 



1/2 



(13) 



Substituting © and (dUJ) in (J7J), and using (|13p gives 



p(u) = - / dq q 2 C(q)w(q) [ 7 + 3(1 - j)P^(q)] 



xe 



[w(q)p/E*] 2 ^P_^ 

du 



or 



du 



V* / dq q 2 C(q)w(q) [7 + 3(1 - j)P*(q)] 

J Qo 



Xe -[w(q) P /E*] 2 ^P 

V 



Integrating this from u = (complete contact, corre- 
sponding to p = oo ) to u gives 



u = v/tt / dq q 2 C(q)w(q) 

J Qo 



/ dp' [ 7 + 3 (l- 7 )J^(g)] 



. (14) 



Let us consider the limiting case of very low squeezing 
pressure. If we introduce x = w(q)p f / E* the last integral 
in (p3J) becomes 

J= dx - [ 7 + 3(1 - j)P 2 (x)} e~ x \ 

Jpw(q)/E* x 

where 

p{x) = ^- e -^ /2 . 

A Jo 

Performing a partial integration gives 



J = 



logx [ 7 + 3(l-7)P 2 (x)]e-^ 



pw(q)/E* 



[ dx \ogx (6(1 - j)P(x)P'(x) 

Jpw(q)/E* V 



+ [ 7 + 3(l- 7 )P 2 (*)] (-2x))e- x \ 
The leading contributions to J as p — > is 

J = "7log (^r) " J™ dx logx (6(1 - 1 )P{x)P\x) 

+ [ 7 + 3(1 - 7 )P 2 (x)] (-2x))e- 2 = - 7 log . 

(15) 



where 



e = exp 



.1^7 

7 



dx logx (q^—-±P{x)P'{x) 



l + ??—tp 2 {x) 

7 



(-2a;)) t 



(16) 



Using (16) we get e = 4.047. Note that e depends on 
Pp(^) for all pressures from small relative contact area 
(P << 1) to complete contact corresponding to P = 1. 
Thus, although the slope of the linear relation between 
u and logp, which hold for very small p, only depend on 
P p (q) for very small p (where the relative contact area 
is proportional to p), the lateral position of the line does 
depend on the whole function P p (q) [or P(x)]. For this 
reason it is important to accurately describe how U e \ de- 
pend on P p (q) for all even if one is only interested in 
the relation between u and p for very small p. The phys- 
ical reason for this is simple: even for arbitrary small 
applied nominal stress p, the stress (at high enough mag- 
nification) in the area of contact will be very high, which 
may result in (nearly) complete contact in the asperity 
contact regions. 

Substituting (15) in (14) gives 



J Qo 



dq q 2 C(q)w(q)>ylog 



eE* 



-u [logp-log(f3E*)}, 



or 



where 



and 



p = f3E*e~ 



a/u 



U 



f3 = eexp 



/ dq q C(q)w(q), 

J QO 



f£dq q 2 C(q)w(q)\og[w(q)} 
r dq q 2 C(q)w(q) 



(17) 



(18) 



(19) 



The relation dSj) between p and u for the special case 
of self-afflne fractal surfaces was studied in Ref. [4] using 
W(q,p) = / yP p (q) in the expression for the elastic energy. 

4. Molecular Dynamics 

We have performed Molecular Dynamics (MD) to study 
the contact area and the interfacial separation from small 
contact to full contact. We are interested in surfaces with 
random roughness with wavelength components in some 
finite range Ao > A > Ai , where Ao is typically similar to 
(but smaller than) the lateral size of the nominal contact 
area. In order to accurately reproduce the contact me- 
chanics between elastic blocks, it is necessary to consider 
solid block which extends at least a distance ~ Ao in the 
direction normal to the nominal contact area. (Note: the 
lower part of the system has been called substrate while 
the upper part is called block.) This leads to an enor- 
mous number of atoms or dynamical variables even for a 
small systems. In order to avoid this trouble we have de- 
veloped a multiscale MD approach. This approach have 




FIG. 3: Surface height profile of a mathematically generated 
self afflne fractal surface (104 x 104 nm 2 square surface area) 
with the root mean square roughness 1 nm. (a) high magni- 
fication (C = 216), (b) low magnification (£ = 4) 

been described in detail in Ref. [l(j and is only summa- 
rized here. The system has lateral dimension L x = N x a 
and L y = N y a, where a is the lattice space of the block. 
Periodic boundary condition has been used in xy plane. 
For the block N x = N y = 400, while the lattice space 
of the substrate b « a/0, where <fi = (1 + y/E)/2 is the 
golden mean, in order to avoid the formation of commen- 
surate structures at the interface. The mass of the block 
atoms is 197 a.m.u. and the a = 2.6 A. The elastic mod- 
ulus and Poisson ratio of the block is E = 77.2 GPa and 
v = 0.42, respectively. 

The atoms at the interface between the block and 
substrate interact with the repulsive potential U(r) = 
e (ro / r) , where r is the distance between a pair of 
atoms. We use r = 3.28 A and e = 74.4 meV. In the 
MD-model calculations there is no unique way of how to 
define the separation u between the solid walls (see Ref. 
[Tot for a discussion of this point). We have used the 
same definition as in Ref. [10[ u = d — d c , where d is 
the separation between the average z coordinate of the 
bottom layer of the block atoms and the average plane of 
the substrate. d c is the critical atom-atom separation we 
use to define contact on atomic scale. Thus, u = cor- 




FIG. 4: 3D images of contacting regions with both substrate 
and block, at low magnification £ = 4 under different load 
(a) ao/E* = 0.0129 and A/A = 0.1089, (b) a /E = 0.26 
and A/A = 0.9527. Note E* = E/(l - is 2 ) is an effective 
modulus. 
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FIG. 5: The probability distribution Pu of surface height h for 
the mathematically generated surface shown in Fig. [3] For a 
square area 104 x 104 nm 2 with lattice constant a = 1.605 A. 
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FIG. 6: Surface roughness power spectra C(q) of the mathe- 
matically generated substrate surface studied in Ref. [19[ and 
below. For a self affine fractal surface with the fractal dimen- 
sion Df = 2.2 and q L = 6.04 x 10 7 m" 1 , q = 1.81 x 10 8 m -1 , 
and qi = 3.9 x 10 10 m -1 , and the root mean square roughness 
amplitude /i rms — 1 nm. 

responds to the separation d c = 4.3615A between planes 
through the center of mass of the interfacial atoms of the 
block and the substrate. 

5. Numerical results: Contact area 

From our molecular dynamics simulation we can calcu- 
late the interfacial stress distribution. In order to obtain 
the contact area we follow the procedure outlined in Ref. 
[Tot and fit the numerical results to the theoretically pre- 
dicted stress distribution 

y ' 2(ttG) i /2 v ; 

(15) 

where G(p, Q depends on the nominal squeezing pressure 
p and the magnification £ (but which is independent of 
cr), and which we choose to get the best agreement with 
numerical data. In Fig. 0we have shown the good agree- 
ment between the numerical pressure distribution and 
the analytical function ([T5]) for £ = 4 under three differ- 
ent nominal pressure. When G is known we can calculate 
the relative contact area using 

— = / daP(aX), (21) 

which is equivalent to using (4) with G replaced by G. 
In Fig. [8] we show the result for A/Aq as a function of 
normalized pressure p/E* . 

In FigE] we have compared the contact area between 
MD simulations and continuum contact mechanics the- 
ory of Persson. Note that the simulation for small load 
predict a contact area which is about ~ 30% larger than 
predicted by the theory. This is slightly larger than what 
has been found in earlier numerical simulations. Thus, 
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FIG. 7: The pressure distribution for £ = 4 for three different 
nominal pressure. The analytical theory has been used to fit 
the numerical pressure distribution as good as possible. Note 
r = G/G (see Ref.[ljj). 



the finite element calculations of Hyun and Robbins[28[ 
and the Green's function molecular dynamics study of 
Campana and Miiser[ll] give a contact area which is 
about ~ 20% larger than that predicted by the Persson 
theory. Similarly, the study of H6nig[29| gives a contact 
area which is about ~ 25% bigger than the analytical one 
for small load. However, none of the computer simula- 
tions can be considered as perfectly converged, and the 
difference between theory and fully converged computer 
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FIG. 8: Contact area ratio A/Ao calculated from (21), as a 
function of normalized pressure p/E* . 
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FIG. 9: Projected contact area comparison between molecular 
dynamics simulations Amd and continuum mechanics theory 
of Persson Ath- At small squeezing pressure, Amd is about 
(30 — 38)% bigger than the Ath- However, the difference 
decreases with increasing pressure. 
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FIG. 10: An elastic block squeezed against a rigid rough sub- 
strate. The (natural) logarithm of the normalized average 
pressure p/E* , as a function of the normalized separation be- 
tween the average plane of the substrate and the average plane 
of the lower surface of the block and the average plane of the 
lower surface of the block denoted by u/h r ms- Adapted from 
0. 
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simulations may be smaller than that indicated by the 
numbers given above. Thus, most of the numerical stud- 
ies reported use rather few grid points within the shortest 
wavelength roughness, which results in an overestimation 
of the contact area [30[. In our simulation for ( = 4 we 
have many atoms within the shortest substrate roughness 
wavelength, but the surface roughness extends over less 
than one decade in length scale. 

Finally, we note that while the pressure distribution we 
obtain for low magnification (( = 4) is in a good agree- 
ment with Persson's theory, for the highest magnification 
(( = 216) this is not the case because only one atom (or 
of order one atom) occurs within the shortest wavelength 
roughness of the substrate (which is roughly given by the 
substrate lattice constant). In the latter case we observe 
that the stress probability distribution for high normal 
stress falls of roughly exponentially rather than like a 
Gaussian. This has also been observed in earlier (non- 



FIG. 11: The relation between the (natural) logarithm of the 
squeezing pressure p (normalized by E*) and the interfacial 
separation u (normalized by the root-mean-square roughness 
amplitude h rms ) for an elastic solid squeezed against a rigid 
surface with the power spectra shown in Fig. [3] In the calcu- 
lation we have used 7 = 0.42. Adapted from (13T | . 

converged) finite element calculations [9[ . It is clear that 
this limiting case cannot be described by the elastic con- 
tinuum model, which is the basis for all analytical contact 
mechanics theories. 

6. Numerical results: Interfacial surface sepa- 
ration 

In Fig. [10] we show the (natural) logarithm of the 
normalized average pressure p/E*, as a function of the 
normalized separation u/h rms between the average plane 
of the substrate and the average plane of the lower surface 
of the block. We show results for the magnification £ = 4 



(open circles) and £ = 216 (solid squares). In this figure 
u = corresponds to the separation 4.3615A between 
the plane through the center of the atoms of the top 
layer of substrate atoms and bottom layer of block atoms. 
Since the atoms interact with a long-range repulsive oc 
r -12 pair potential, it is possible to squeeze the surfaces 
closer to each other than what corresponds to u = 0. 
This explains why simulation data points occur also for 
u < 0. The theory described in Sec. 3 assumes a contact 
interaction potential so that u > 0, and can therefore 
not be compared with the MD simulations for very small 
(and negative) u. 

In Fig. [TT] we compare the MD results from Fig. [10] 
(solid squares) with the theoretical prediction calculated 
from (p3j) using the same surface roughness power spectra 
(and other parameters) as in the MD-calculation. The 
theory is in good agreement with the numerical data for 
0.2 < u/h rms < 2. For u/h rms < 0.2 the two curves dif- 
fer because of the reason discussed above, i.e., the "soft" 
potential used in the MD simulation allows the block and 
substrate atoms to approach each other beyond u = 0, 
while in the analytical theory a contact potential is as- 
sumed where the repulsive potential is infinite for u < 
and zero for u > 0. The difference between the theory 
and the MD results for u/h rms > 2 is due to a finite 
size effect. That is, since the MD calculations use a very 
small system, the highest asperities are only ~ 3h rms 
above the average plane (see the height distribution in 
Fig. [5]), and for large u very few contact spots will occur, 
and in particular for u > 3h rrns no contact occurs and 
p must vanish. In the analytical theory, the system is 
assumed to be infinite large. So that even for a Gaussian 
distribution of asperity height, there will always be (in- 
finitely) many infinitely high asperities and contact will 
occur at arbitrary large separation u. The asymptotic re- 
lation u oc \ogp will hold for arbitrarily large u (or small 
squeezing pressures p) . 

7. Contact mechanics for a measured surface 

Pei et al[30] have performed a Finite Element Method 
(FEM) computer simulation of the contact mechan- 
ics for a polymer surface, using the measured surface 
topography [24J as input, squeezed against a flat surface. 
Here we would like to compare the FEM results with the 
analytical results of Persson. 

In Fig. [12] we show the surface height profile of a poly- 
mer surface studied in Ref. [24] with the root mean 
square roughness 14 nm. The probability distribution 
of surface height P^, for two different square 10 x 10 /xm 2 
surface areas on the polymer film, has been shown in 
Fig. [HI These two surface areas were used in the FEM 
calculation by Pei et al as the upper and lower surface. 
It is remarkable that in spite of the fact that the height 
distributions are not perfect Gaussian (in particular one 
surface exhibits a "bump" in the height distribution on 
the outer side of the height profile), Pei et al obtained a 




FIG. 12: Surface height profile of a polymer surface, with the 
root mean square roughness 14 nm, measured over a 10 x 
10 /im 2 square surface area. 
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FIG. 13: The probability distribution Pu of surface height h 
for the polymer film in Fig. 1121 The two curves correspond to 
height profile data measured at two different square surface 
areas, each 10 x 10 /xm 2 . 

nearly perfect linear relation between \ogp and u. This 
result indicates that even in the present case the area 
of real contact is proportional to the load (see below) 
and the statistical properties of the contact regions does 
not change with the load for small load. The surface 
roughness power spectra of the two surfaces are shown in 
Fig. [EJ and are almost identical in spite of the difference 
which occurs in the height distribution P^. 

Fig. [T5l shows the relation between the logarithm (with 
10 as basis) of the squeezing pressure p (normalized by 
E*) and the interfacial separation u (in nm) for an elas- 
tic solid squeezed against a rigid surface with the power 
spectra given by the sum of the two power spectra shown 
in Fig. [H] The line a is the result of the FEM calcu- 
lation of Pei et al, and shows that for large separation 
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FIG. 14: Surface roughness power spectra C(q) of the polymer 
film in Fig. [12] The two curves correspond to height profile 
data measured at two different square surface areas, each 10 x 
10 /im 2 . The height probability distribution for the same 
surface areas is shown in Fig. [13] 
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FIG. 15: The relation between the logarithm (with 10 as ba- 
sis) of the squeezing pressure p (normalized by E*) and the 
interfacial separation u (in nm) for an elastic solid squeezed 
against a rigid surface with the power spectra given by the 
sum of the two power spectra shown in Fig. [14] The line a is 
the result of the Finite Element calculation of Pei et al, while 
the line b is the prediction of the theory. In the calculation 
we have used 7 = 0.38. 



p oc exp(— u/jhrms) with 7 « 0.4, which is consistent 
with our analytical results (see line b in Fig. [T5|) . 

Fig. [16] shows the relation between the relative con- 
tact area, A/Aq, and the normalized squeezing pressure, 
p/E* , for an elastic solid squeezed against a rigid surface 
with the power spectra, given by the sum of two power 
spectra shown in Fig. [14] The line a is the result of the 
FEM calculation of Pei et al, while the line b is the pre- 
diction of the theory where A/Aq is scaled by a factor of 
1.29. 

The good agreement, between the FEM calculations 




A/Ao 

FIG. 16: The relation between the relative contact area, 
A/Ao, and the normalized squeezing pressure, p/E*, for an 
elastic solid squeezed against a rigid surface with the power 
spectra given by the sum of the two power spectra shown in 
Fig. 1141 The line a is the result of the Finite Element cal- 
culation of Pei et al, while the line b is the prediction of the 
theory where A/Aq is scaled by a factor of 1.29. 



and the analytical theory found above, indicates that the 
contact mechanics results are robust and not very sensi- 
tive to many details such as the assumption of perfectly 
random surfaces, which is unlikely to be exactly obeyed 
for the polymer surfaces, the surface topography of which 
was used in the FEM calculations. (In order to address to 
what extent a measured surface is randomly rough, one 
would need to calculate higher order correlation func- 
tions, and show that odd order functions in the height 
coordinate vanish (or are very small), and that even or- 
der height correlation functions can be decomposed into 
a sum of pair correlation functions. We are aware of no 
such study for "real" measured surface profiles.) 

Finally we note that the observation of an effective 
exponential repulsion has important implications for tri- 
bology, colloid science, powder technology, and materials 
science [24|. For example, the density or volume of granu- 
lar materials has long been known to have a logarithmic 
dependence on the externally applied isotropic pressure 
or stress, as found, for example, in the compression stage 
during processing of ceramic materials [3l[. Recent work 
on the confinement of nanoparticles has also indicated an 
exponential force upon compression [32], suggesting that 
this relationship could be prevalent among quite different 
types of heterogeneous surfaces. 

8. Variation of the average surface separation 
u(() with the magnification £ 

The theory presented above can be easily generalized 
in various ways. Thus, it is possible to include the ad- 
hesional interaction [3, [33[. In this case the work done 
by the external pressure p will be the sum of the stored 



magnification £ 

hard solid: 




FIG. 17: An asperity contact region observed at the magnifi- 
cation f . It appears that complete contact occur in the asper- 
ity contact region, but upon increasing the magnification it is 
observed that the solids are separated by the average distance 

a(C). 




FIG. 18: The dependence on the magnification, of the the av- 
erage distance u(() (in units of the root-mean-square rough- 
ness amplitude forms (1) of the whole surface) between the sur- 
faces in an asperity contact region observed at the magnifica- 
tion £. For the polymer surface with the power spectra shown 
in Fig. [14] We have assumed the effective elastic modulus 
E* = 2 GPa and the nominal squeezing pressure a: po = 1, 
b: 10, c: 100 and d: 200 MPa. 



(asperity induced) elastic energy plus the (negative) ad- 
hesional energy, i.e., the right hand side of ([6]) will now 
be U G \ + E/ad- The theory can also be applied to study 
how the spacing u(Q depends on the magnification. Here 
u(() is the (average) spacing between the solids in the ap- 
parent contact areas observed at the magnification £. 
Fig. [T71 shows an asperity contact region at the magni- 
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FIG. 19: The dependence on the magnification, of the rel- 
ative contact area A(£)/Ao, the root-mean-square roughness 
amplitude forms (C) ( m units of the root-mean-square rough- 
ness amplitude forms (1) of the whole surface) and the average 
distance u(Q between the surfaces in an an asperity contact 
region observed at the magnification All quantities are 
shown on a logarithmic (with 10 as basis) scale. For the poly- 
mer surface with the power spectra shown in Fig. 1141 We 
have assumed the effective elastic modulus E* = 2 GPa and 
the nominal squeezing pressure po = 10 MPa. 



fication £. It appears that complete contact occurs in the 
asperity contact region, but upon increasing the magni- 
fication it is observed that the solids are separated by 
the average distance u(Q. The information about how 
u(() depends on the magnification is crucial for many 
important applications, e.g., sealing (see below). 

We study the contact between the solids at increas- 
ing magnification. In an apparent contact area observed 
at the magnification £, the substrate has the root mean 
square roughness amplitude 

rqi 

h 2 rms (C) = 2tt / dq qC(q) 

The separation between the surfaces u(Q is given by (14) 
but with the lower integration limit given by qo( instead 
of qo and the lower pressure integration limit given by 
p(0=PoA /A(0: 



(0 = vW dq q 2 C(q)w(q) 
JCqo 



x 



/ dp' i [ 7 + 3(1 - 7 )p 2 (<?y , o] e- [wiqX)p /E1 

J via p 



MO 
where 



(23) 



w(q,Q = (* J" dq' q ,3 C(q'^j (24) 

and where P(q,p f ,() is given by (12) with s — 
w(q,Q/E*. When we study the apparent contact area 



at increasing magnification, the contact pressure p(Q will 
increase and the surface roughness amplitude h rms (Q will 
decrease. Thus, the average separation u(Q, between the 
surfaces in the (apparent) contact regions observed at the 
magnification £, will decrease with increasing magnifica- 
tion. 

In Fig. IT8I we show, for the polymer surface with the 
power spectra shown in Fig. [I4j how the logarithm of the 
average distance u(() (in units of the root-mean-square 
roughness amplitude h rms (l) of the whole surface), de- 
pends on the magnification £. We have assumed the ef- 
fective elastic modulus E* = 2 GPa and the nominal 
squeezing pressures a: pq = 1, b: 10, c: 100 and d: 
200 MPa. Note that for ( > 2 the two lowest squeez- 
ing pressures give virtually identical separation between 
the surfaces in the (apparent) asperity contact regions, 
in spite of the fact that the pressure for curve b is 10 
times higher than for curve a. This just reflects the fact 
stated earlier that for low squeezing pressure the area 
of (apparent) contact A varies linearly with the squeez- 
ing force pAo, and the interfacial stress distribution, and 
the size-distribution of contact spots, and the interfa- 
cial separation u((), are independent of the squeezing 
pressure H, [5^ . That is, with increasing p, existing con- 
tact areas grow and new contact areas form in such a way 
that in the thermodynamic limit (infinite-sized system) 
the quantities referred to above remain unchanged. 

In Fig. [19] we again show how the average distance 
u(() depends on the magnification (. We also show the in- 
dependence of the relative contact area A(Q/Aq, and the 
root-mean-square roughness amplitude forms (C) ( m units 
of the root-mean-square roughness amplitude h rms (1) of 
the whole surface). All quantities are shown on a log- 
arithmic (with 10 as basis) scale, and the results are 
for the polymer surface with the power spectra shown 
in Fig. [TH We have assumed the effective elastic mod- 
ulus E* = 2 GPa and the nominal squeezing pressure 
po = 10 MPa. Note that for ( < 2, u(()/h rms (l) « 1.23 
i.e., at the lowest magnification the upper surface is "rid- 
ing" on top of the largest substrate asperities. The sepa- 
ration between the solids in the asperity contact regions 
rapidly drop with increasing magnification (correspond- 
ing to smaller and smaller asperity contact regions), and 
already at the magnification ( = 10 the separation is only 
~ 3% of the average separation which occurs at the low- 
est magnification. This is, of course, mainly due to the 
strong increase in the local pressure (as manifested in the 
decreased contact area, A(10)/Aq « 0.1) in the asperity 
contact regions at high magnification but also due to the 
reduction in the effective roughness detected over short 
distances (the largest contribution to h rms (1) comes from 
the longest wavelength roughness components). 

Let ui(() be the (average) height separating the sur- 
faces which appear to come into contact when the mag- 
nification decreases from ( to ( — A£, where A£ is a 
small (infinitesimal) change in the magnification. The 
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FIG. 20: The interfacial separation u((^) and i^i(C) as a func- 
tion of the magnification £ for the surface shown in Fig. [3] 
and for the squeezing pressure 1.3 GPa. 
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FIG. 21: The probability distributions P u and P u are defined 
in the text. For the surface shown in Fig. [3] and for the 
squeezing pressure 1.3 GPa. The distribution P u has a delta 
function at u = with the weight A((i)/Aq « 0.047. 



empty volume between the surfaces which appears to 
be in contact at the magnification £ — A£ is given by 
u(( — A()A(( — A£). But this volume must be the sum 
of the volume u(QA(Q between the surfaces which ap- 
pears to be in contact at the magnification plus the 
additional volume Ui(()[A(( — AQ — A(Q] which results 
from the increase in the apparent contact area as the 
magnification decrease from £ to £ — A£: 

m(C-AC)A(C-AC) = u(C)A(0+ Ul (C)[A(C-AO-A(C)] 



or 



«i(C) = «(C) + « , (CM(C)M'(C) 



(25) 



In Fig. [20] we show the interfacial separation u(Q and 
u\ (£) as a function of the magnification £ for the surface 
shown in Fig. [3] and for the squeezing pressure 1.3 GPa. 



Note that for small magnification u\ increases rapidly. 
This can be understood as follows: The theory is for an 
infinite system. For an infinite system, even for a Gaus- 
sian distribution of surface height, there will always be 
some infinitely high asperities and some infinitely deep 
valleys. Thus during contact there will always be some 
regions where the surface separation is arbitrarily high. 
Of course, the fraction of the surface where u\ is large is 
extremely small. The reason is that the variation of the 
(apparent) contact area with the magnification is negli- 
gible in the region where u\ start to grow fast (see Fig. 
[20]) . Thus, the strong increase in u\ for small magnifica- 
tion is of no practical importance-it is a purely academic 
result. 

The quantity ui(() is very important in the context of 
the leakage through rubber sealing: Let us study the in- 
terface between the rubber and the substrate as the mag- 
nification increases. At low magnification it appears as if 
the solids makes perfect contact at the interface. As we 
increase the magnification non-contact area becomes vis- 
ible. At large enough magnification, say ( = £*, the non- 
contact area will percolate [19]. A first rough estimate of 
the gas (or fluid) leakage is obtained by assuming that 
the gas flow through a pipe or pore with width and length 
A « L/(* (where L is the linear size of the sealing con- 
tact area) and height u±(C*), and that the whole pressure 
drop Ap = pi—po (where p\ and po is the pressure to the 
left and right of the sealing) occurs over the pore. Thus 
for an incompressible fluid, the mass-flow per unit time 
through the interfacial pore will be Q « pu\(C ) *)Ap/12r] 
(where rj is the viscosity). We will analyze this problem 
in greater detail elsewhere [34|. 

Finally, let us consider the distribution of interfacial 
separations [13] 

P u = (S(u - u(x))> (26) 

where (..) stands for ensemble average (which often is 
equivalent to average over the surface area) . We also de- 
fine another distribution P(u) of interfacial separations 
which differ from (26) by using instead of u(x) another 
function which involves some average over the spatial co- 
ordinate and defined as follows: The probability to find 
the surface separation u < ^i(C) is 

*(«) = A(C)/A = P(() 

where £ = ((u) is the solution to u = u\{Q [note: u\(Q 
is a monotonically decreasing function of £ so there exists 
a unique solution £ = ((u) to u = i/i(C)]- The probability 
distribution 

- _ dV(u) _ P'(Q , 

du ~ «i(C) C=C(U) 

We can also write 

P u = -Jd(^-5[u- Ul (0} (27) 



It is convenient to change integration variable to /i de- 
fined by ( = exp(— n) and consider A and Wi as function 
of ji. This gives 

Pu = - [ d»^-5[u- Ul ([i)} (28) 

In Fig. [21] we show P u and P u for the surface shown in 
Fig. [H and for the squeezing pressure 1.3 GPa. The re- 
sult for P u was obtained from the MD-simulations, while 
P u was obtained from the analytical theory presented 
above. As expected the distribution P u is more narrow 
than P u [since it involves u\ which is already an aver- 
age of u(x)], but it is easy to show that average of u 
is the same for both P and P and equal to the average 
separation between the surfaces. Thus 

f duuP u = - f d( ^-MO 

Substituting (25) in this equation gives 

JduuP u = -±- £ d( ±[A(()u(0] = «(Co) 

where we have used that ^(Ci) = and that A((o) = Aq 
(note: ("o is the lowest magnification usually chosen as 
unity). Since u((o) equals the average separation between 
the surfaces we have proved our statement. Note that P u 
has a delta function at u = with the weight A((i)/Aq. 
Using this fact it is easy to show that distribution P u has 
the zero order moment equal to unity: 

J duP u = l 

Thus, the zero order and first order moments of P u and 
P u are the same, but higher order moments will differ. 

9. Elastic energy and adhesion 

In Sec. 3 we have shown that from the dependence 
of the surface separation u on the squeezing pressure p, 
deduced from, e.g., experiments or from MD or FEM 
calculations, one can obtain the elastic energy U e \ stored 
in the asperity contact regions. Thus, in Sec. 3 we have 
presented an expression for U e \ which results in nearly 
the same relation between u and p as observed in the 
numerical MD and FEM studies. This is an important 
result as U e \ is relevant for many important applications, 
e.g., adhesion between elastic solids with randomly rough 

surfaces d El & 

We consider the adhesive contact between two elas- 
tic solids with randomly rough surfaces. Assume that 
the surface roughness power spectra has a long-distance 
roll-off wavelength Ao (corresponding to the roll-off wave 
vector qo = 2tt/\o) which is much smaller than the di- 
ameter of the nominal contact area. In this case we can 
take into account the influence of the surface roughness 
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FIG. 22: The effective interfacial binding energy and the area 
of real contact as a function of the normalized root-mean- 
square roughness /irms/formso- The case h rms — /irmsO COr- 
responds to the power spectra obtained from the measured 
height profile for the surface 7 (with h Ym so ~ 0.2 fim) studied 
in Ref. [37j] (see also Fig. [23]) . In the calculation we have used 
the measured (low- frequency) elastic modulus (E « 5 MPa) 
and the measured (for flat surfaces) interfacial binding energy 
(per unit area) (A7 ^0.1 J/m 2 ). In (a) we show the calcu- 
lated 7eff (1) (red curve) and the atomic contact area A(£i)/Aq 
observed at the highest magnification. The green curve gives 
7eff(l) under the assumption that complete contact occurs at 
the interface. When the system is in thermal equilibrium it 
will be in the state where 7 e ff is maximal. In (b) we show the 
thermal equilibrium interfacial binding energy and the corre- 
sponding contact area. Note that the system flips abruptly 
from the complete contact to partial (about 50%) contact at 

hrms ~ 0.9/lrmsO- 

on the (adhesive) contact mechanics by using an effective 
interfacial binding energy 

7eff(CM*(C) = A 7 A*(Cl)-t/el(C), 

where U e \(() is the elastic energy stored in the asperity 
contact region as a result of the asperities which can- 
not be observed at the magnification ( (i.e., due to the 
surface roughness with wavelength shorter than \o/(). 
A*(C) is the contact area when the surface is studied at 
the magnification which in general is larger than the 
projected contact area A(Q. The interfacial binding en- 
ergy per unit surface area for the contact between two 
perfectly flat surfaces of the two solids is denoted by A7, 




root-mean-square roughness (jam) 

FIG. 23: The normalized effective interfacial binding energy 
7e/ / / A7 as a function of the root-mean-square roughness for 
7 differently prepared surfaces. Blue curve: theory, assuming 
perfect (atomic) contact at the interface; green curve: a; and 
red curve: b. The experimental data for the pull-off velocity 
0.2 and 2 /rni/s. From ref. [13 



and A*((i) is the contact area observed at the highest 
(atomic) magnification £i- Using (9) we can write [3] 

7..K) = - *-5^j£ <"> M*p)CW 

The macroscopic effective interfacial binding energy 
7 e ff(l) determines the macroscopic contact mechanics 
and the pull-off force. For example, for a rubber ball 
(radius R) in contact with a nominal flat substrate, the 
pull-off force is given by the JKR formula F pu ii-off = 
(37r/2)i?7 e ff. Thus, the surface roughness enters in the 
expression for (1), but the surfaces can be considered 
as perfectly smooth when solving the macroscopic con- 
tact mechanics problem. 

In Ref. [3] one of us studied 7 e ff(C) using the approxi- 
mation W(q,p) = P p (q) when calculating the elastic en- 
ergy. Here we will present results using the improved 
expression for U e \ with W(q,p) given by (10). In Fig. [22] 
we show the effective binding energy and the area of real 
contact as a function of the normalized root-mean-square 
roughness h rms /h rms o. The case forms — formsO corre- 
sponds to the power spectra obtained from the measured 
height profile for the surface 7 (with fo rms o ~ 0.2 fim) 
studied in Ref. [13 (see also Fig. [23]). The power spectra 
used in Fig. [22] have been obtained by scaling the power 
spectrum of surface 7 with the factor ( h rms /h rms o) 2 . In 
the calculation of 7 e ff we have used the measured (low- 
frequency) elastic modulus (E w 5 MPa) and the mea- 
sured (for flat surfaces) interfacial binding energy (per 
unit area) (A7 w 0.1 J/m 2 ). In Fig. [22](a) we show 
the calculated 7 e ff(l) (red curve) and the atomic con- 
tact area A((i)/Aq observed at the highest magnifica- 
tion. The green curve gives 7 e ff (1) under the assumption 



that complete contact occurs at the interface. When the 
system is in thermal equilibrium it will be in the state 
where the interfacial binding energy 7 e fj is maximal. In 
Fig. [22] (b) we show the thermal equilibrium interfacial 
binding energy and the corresponding contact area. Note 
that the system flips abruptly from the complete contact 
to partial (about 50%) contact at h rms « 0.9/i rms o. Very 
similar abrupt transitions have been found in exact solu- 
tion for a cosines roughness profile, see e.g., Ref. [38]. We 
also note that while Fig. l22lf b) shows the minimum free 
energy state, hysteresis may occur in practical situations 
(see Ref. [38]). Finally, note that the contact area A((i) 
is finite also when 7 e ff(l) = 0, i.e., adhesion increases 
the contact area even if no adhesion can be detected in 
a pull-off experiment. Since it is the area of real contact 
which determines the sliding friction force, the adhesional 
interaction may increase the sliding friction even if no ad- 
hesion can be detected in a pull-off experiment. We plan 
to study the adhesive contact between randomly rough 
surfaces using MD simulations to compare with the pre- 
dictions of the theory described above. 

In an earlier publication one of us has studied the ef- 
fective binding energy as a function of the root-mean- 
square roughness for 7 differently prepared surfaces [37j 
(see Fig. [23]) . The theory curve (blue curve) in Fig. [23] is 
based on the assumption of perfect (atomic) contact at 
the interface. The curves a and b are experimental re- 
sults for the pull-off velocity 0.2 and 2 /im/s. The theory 
data point for 7 e ff(l)/A7 for the highest roughness value 
(/i rms = /irmso ~ 0.2 /im) was actually negative (about 
—0.1, see Fig. [22]) . indicating that the complete con- 
tact state cannot be the ground state, and in Fig. ([23]) 
we therefore gave the binding energy 7 e ff(l) = 0, of the 
fully detached state. However, in accordance with the ex- 
perimental data points, we find that with the improved 
elastic energy used above, the partly contact state has a 
lower free energy (a larger binding energy) than the fully 
contact state, i.e., 7 e ff(l) is positive for the surface 7, 
albeit still somewhat smaller than observed experimen- 
tally (compare Fig. [22] for h rms /h rms o = 1 with Fig. [23] 
for h rms0 = 0.2 /im). 

10. Summary and Conclusion 

We have used our recently developed multiscale molec- 
ular dynamics approach [1(| to study real contact area 
and interfacial separation from small contact to full con- 
tact. The real contact area rises linearly with the load 
for small load. Here we have found that, at low magnifi- 
cation where the atomistic nature of the solids becomes 
unimportant, the MD results match very well with Pers- 
son's theory, especially when the contact approaches full 
contact. 

The interfacial separation as a function of squeezing 
pressure has been derived theoretically. For non-adhesive 
interaction and small applied pressure, p oc exp(— u/uo) 
is in good agreement with experimental observation. 



This relation has been tested with MD simulations and 
they match quite well with each other. The present 
results may be of great importance for soft solids, e.g. 
rubber-like material, or very smooth surfaces. 
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